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We probe the possibility that Centaurus A (Cen A) is a point source of ultra-high 
energy cosmic rays (UHECR) observed by PAO, through the statistical analysis of 
the arrival direction distribution. For this purpose, we set up the Cen A domi- 
nance model for the UHECR sources, in which Cen A contributes the fraction fc 
of the whole UHECR with energy above 5.5 x 10 19 eV and the isotropic background 
contributes the remaining 1 — fc fraction. The effect of the intergalactic magnetic 
fields on the bending of the trajectory of Cen A originated UHECR is parameter- 
ized by the gaussian smearing angle 9 S . For the statistical analysis, we adopted the 
correlational angular distance distribution (CADD) for the reduction of the arrival 
direction distribution, and the Kuiper test to compare the observed and the expected 
CADDs. Using CADD, we identify the excess of UHECR in the Cen A direction 
and fit the CADD of the observed PAO data by varying two parameters fc and 9 S 
of the Cen A dominance model. The best-fit parameter values are fc ~ 0.1 (The 
corresponding Cen A fraction observed at PAO is /c,pao ~ 0.15, that is, about 10 
out of 69 UHECR.) and 9 S = 5° with the maximum probability P max = 0.29. Con- 
sidering the uncertainty concerning the assumption of isotropic background in the 
Cen A dominance model, we extend the viable parameter ranges to the 2a band, 
0-09 < /c,pao < 0.25 and 0° < S < 20°. This result supports the existence of a 
point source extended by the intergalactic magnetic fields in the direction of Cen A. 
If Cen A is actually the source responsible for the observed excess of UHECR, the 
average deflection angle of the excess UHECR implies the order of 10 nG intergalactic 
magnetic field in the vicinity of Cen A. 
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I. INTRODUCTION 



The origin of the ultra-high energy cosmic rays (UHECR) is a long-standing puzzle (if. 
The recent confirmation of the Greisen-Zatsepin-Kuzmin (GZK) suppression in the cosmic 
ray energy spectrum [2|-0] implies that UHECR with energies above the GZK cutoff, Eq ZK ~ 
4 x 10 19 eV, mostly come from relatively close extragalactic sources within the GZK radius 
fGZK ~ 100 Mpc. Furthermore, UHECR with these energies are expected not to be strongly 
affected by the galactic or extragalactic magnetic field, so that their arrival directions keep 
some correlation with the source distribution and can be used to trace the sources of UHECR. 
Recently Pierre Auger Observatory (PAO) released the updated UHECR data with energy 
E > 5.5 x 10 19 eV [5[. These data can be used for tracing the distribution of UHECR 
sources through the statistical analysis of the arrival direction distribution. Beginning with 
the report of the correlation between UHECR and active galactic nuclei (AGN) by the 
PAO collaboration 0-0], several attempts was made to identify the UHECR sources. Most 
attempts aimed to test the plausibility of a certain kind of high energy astrophysical objects 
whose number ranges from a few to a few hundred [8-3|. In this paper, we aim to test the 
other extreme possibility that a single dominating source is responsible for the large part of 
observed UHECR. 

The most mentioned candidates for a strong UHECR source are Centaurus A (Cen A) and 
Messier 87 in Virgo cluster, which are active galaxies very close to us. Many people suggested 
that Cen A might be the source of UHECR EIS 

Cen A is located in the southern sky, 
near the center of the exposure region of the PAO experiment. Therefore, the PAO data 
provide the best chance for checking whether Cen A is a strong source of UHECR. Actually, 
the correlation of UHECR with AGN was strengthened by many UHECR observed around 
Cen A. On the contrary, the number of UHECR around M87 is smaller than the expected 
one considering its distance, weakening the claimed AGN correlation. Our purpose is to 
quantify the Cen A contribution in the observed UHECR by PAO and thus try to establish 
the existence of a point source in the Cen A direction. 

For a single dominating source to be able to explain the large portion of UHECR data, the 
intergalactic magnetic fields must play a significant role, spreading UHECR over the large 
region of the sky around the source. Unfortunately, our knowledge about the intergalactic 
magnetic fields is rather poor yet. Modeling the intergalactic magnetic fields and testing it 
has its own uncertainty. For this reason, we choose a different and simpler strategy. We 
parameterize the effect of intergalactic magnetic fields on the deflection of UHECR trajectory 
by the gaussian spreading of UHECR arrival directions around the sources. In addition, we 
also need to consider the contribution from other sources to explain the whole set of observed 
UHECR arrival directions. For this purpose, we introduce one more parameter measuring 
the fraction of Cen A contribution to observed UHECR. Then we search for the values of two 
parameters with which the observed UHECR arrival directions can be plausibly explained 
through the statistical comparison of the arrival direction distributions. 

Statistical c omp arison of the arrival direction distributions can be done in many different 



ways. In Ref. [23|, |24j, we developed the new statistical test methods, whose basic idea is 
that the two-dimensional distribution of arrival directions is reduced to the one-dimensional 
probability distributions, which can be compared by using the well-known Kolmogorov- 
Smirnov test or its variants. We proposed a few reduced one-dimensional distributions 
suitable for the test of correlation between the UHECR arrival directions and the point 
sources. Among them, we adopt the correlational angular distance distribution method and 
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the flux-exposure value distribution method. These methods will be briefly described in 
Sec. EH 



This paper is organized as follows. In section [TTJ we explain the Cen A dominance model 
for the UHECR sources and the details needed for the Monte-Carlo simulations of UHECR 
arrival directions. In section IHfl we briefly introduce our statistical methods for comparing 
two arrival direction distributions. In Section IIVt the results of our analysis are presented. 
We give a few discussions on the results and conclude in section [Vj 

II. THE SINGLE DOMINATING SOURCE MODEL FOR UHECR 

We examine the plausibility of the idea that the Cen A is the dominant sources of UHECR 
through the statistical analysis of the arrival direction distribution of UHECR. For more 
definite interpretation of our analysis, we need to solidify the UHECR source model with 
the Cen A dominance and the methods adopted for statistical analysis. In this section, we 
describe the details of the Cen A dominance model for UHECR sources. 

Cen A is located at a = 201.37°, 5 = —43.02° in the equatorial coordinates and the 
distance is estimated to be 3 — 5 Mpc {25[. We consider Cen A as a smeared point sources of 
UHECR. This is mainly to incorporate the fact that the trajectories of UHECR can be bent 
by intervening magnetic fields. We may model the intervening galactic and extragalactic 
magnetic fields between Cen A and the earth. But, it would involve large arbitrariness due 
to our lack of knowledge on extragalactic magnetic fields. Instead of such detailed modeling 
of magnetic fields, we simply assume that Cen A has a gaussian flux distribution on the sky 
with a certain angular width 9 S , so that the effect of magnetic fields is measured through the 
parameter 9 S , called the smearing angle. Then the UHECR flux at a direction r contributed 
by Cen A can be written as 



where Fqa is the averaged flux of Cen A, 9(i) = cos _1 (f ■ tca) is the angle between the 
direction r and Cen A, and N(6 S ) = (1/Att) J dQexp[—(9(i)/9 s ) 2 ] is the normalization of 
smearing function. For small 9 S , N(6 S ) « 9 2 /4 and for large 9 S , N(6 S ) « 1. For the small 
smearing angle 9 S , the average deflection angle is (9) = f 9e~^ e ^ 3 ^ dVt/ f e~^ 9 ^ 9s ^ dD, ~ 9 S , 



Though Cen A can be a dominant source of UHECR, it is very unlikely that Cen A is 
the only source of UHECR. We need to consider the contribution from other distributed 
sources. We consider, for the sake of simplicity, that a certain fraction of UHECR are 
originated from Cen A, while the remaining fraction of them are from the isotropically 
distributed background contributions. Then, the expected flux at a give arrival direction r 
is given by the sum of two contributions, 



F CA (f) = F CA 



exp[- (9 J (r)/9 s ) 2 } 
N(9 S ) 



(1) 



and A9 = y/(9 2 } - (9) 2 « 9 s /2. 



F{r) = F CA (f) + F lso . 



(2) 



Now we define the Cen A fraction fo to be 



fc = = 



Fqa + Fiso 



(3) 



i 



The UHECR flux can be written as 

F{i) = fc j ^-m/e,f ] +(i _ fc) - (4) 

where F = F c + -Piso- Out of three parameters F c , F ISO , and 9 S , the Cen A fraction f c 
and the smearing angle 9 S are treated as the free parameters of the model, while the average 
flux F is fixed by the total number of UHECR events. 

To do the simulation for the observed arrival directions of UHECR, we also need to take 
into account the efficiency of the detector as a function of the arrival direction. It depends 
on the location and the characteristics of the detector array. Here we consider only the 
geometric efficiency which is determined by the location and the zenith angle cut of the 
detector array. Then the exposure function h(r) depends only on the declination 5, 

h(S) — — [sin a m cos A cos 5 + a m sin A sin 5] , (5) 

71 

where A is the latitude of the detector array, 8 m is the zenith angle cut, and 



0, for f > 1, 

Oim= { TT; fc> r £ < — 1; With £ 



cos 6 m — sin A sin 5 



cos" 1 ^ otherwise ' cos A cos 5 

The expected flux at the detector array is proportional to F(i)h(i). We also note that the 
Cen A fraction fc is the fraction of Cen A contribution over the whole sky. It is in general 
different from the fraction of Cen A contribution within the sky covered by a given detector 
array, because the latter is masked by the exposure function. We denote the latter, e.g. for 
PAO, by /c,pao to distinguish it from the former. 



III. STATISTICAL COMPARISON OF TWO ARRIVAL DIRECTION 

DISTRIBUTIONS 

We now turn to the statistical methods to measure the plausibility of the Cen A domi- 



nance model. In Refs. [23|, [24j, we developed the simple comparison method for the UHECR 
arrival direction distributions, where the two-dimensional UHECR arrival direction distri- 
butions on the sphere is reduced to one-dimensional probability distributions of some sort, 
so that they can be compared by using the standard Kolmogorov-Smirnov (KS) test or its 
variants. In this paper, we adopt the reduction methods called the correlational angular dis- 
tance distribution (CADD) and the flux-exposure value distribution (FEVD). For a detailed 



explanation on these methods, see Refs. |23|, |24j . Here we present briefly the basic ideas of 



these distributions and how to calculate the probability measuring the plausibility. 

• Correlational Angular Distance Distribution This is the distribution of the angular 
distances of all pairs UHECR arrival directions and the point source directions: 

CADD : {cos 9 if = r< ■ | i = 1, . . . , N; j = 1, . . . , M) , (6) 



where are the UHECR arrival directions, ?'• are the point source directions, and N and 
M are their total numbers, respectively. 
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• Flux Exposure Value Distribution At a given arrival direction, the expected flux value is 
the product of the UHECR flux expected from the UHECR source model and the exposure 
function of the detector at that direction. FEVD is the distribution of expected flux values 
at UHECR arrival directions: 

FEVD : {Fi = F(r,)/i(r 4 ) | i = 1, . . . , N} , (7) 

where r$ are the UHECR arrival directions, N is the total numbers of UHECR, F(ii) and 
h(ii) are the UHECR flux and the exposure function, respectively. 

From the given point source set and the UHECR arrival direction data set, we can 
get CADD and FEVD. Then, we can apply KS test or its variants such as Kuiper test 
and Anderson-Darling test to compare two CADDs or two FEVDs, one from the observed 
UHECR data set and the other from the expected (simulated) UHECR data set from the 
model under consideration. In this analysis, we use Kuiper (KP) test because its sensitivity 
is found to be most appropriate for our purpose and the probability function of its statistic is 
available in analytic form. The KP test is based on the cumulative probability distribution 
(CPD), Sn(x) = f x p(x')dx' and the KP statistic is the sum of maximum difference above 
and below two CPDs, 

D KP = max [S Nl (x) - S N2 (x)] + max [S N2 (x) - S Nl (x)] . (8) 

X X 

From the KP statistic -Drp, the probability that CADD/FEVD of the observed data is 
obtained from the model under consideration can be estimated using the Monte-Carlo sim- 
ulations in general. For the KP statistic -Dkp ; when the data in the distribution are all 
independently sampled, the following approximation formula is available: 

P(D KP \N e ) = g K p([v / A 7 e + 0.155 + 0.24/v/iVepKp), (9) 

where Q K p(A) = 2 Y™ =x {Af\ 2 - V^e' 2 ^ and N e = N 1 N 2 /(N 1 + N 2 ) is the effective number 
of data. For FEVD, the number of data in the distribution is same as the number of UHECR 
data. For CADD, the number of data in the distribution is the number of UHECR data 
times the number of point sources. Therefore, for a single source case, the number of data 
in the distribution CADD and FEVD is same as the number of UHECR data. This means 
that the data in CADD and FEVD are all independent, and we can use the formula (jH]) for 
both CADD and FEVD for a single source model. Now, N\ = No, the number of observed 
UHECR data and N 2 — N$, the number of mock UHECR data. We can make the expected 
distribution more accurate by increasing the number of mock data N 2 . In the limit N 2 — > oo, 
the effective number of data is simply N e = No- 

IV. THE RESULTS OF STATISTICAL TESTS 

For the observed UHECR data set, we use the UHECR data released by PAO in 2010 (j. 
The released data set contains 69 UHECR with energy higher than 5.5 x 10 19 eV. The PAO 
site has the latitude A = —35.20° and the zenith angle cut of the released data is 8 m = 60°. 
The arrival direction distribution of the released PAO data is shown in Fig. [TJ The locations 
of Cen A and M87 are also marked for reference. For the correlation test between the 
UHECR arrival directions and the astrophysical objects, the choice of the energy cut can be 
crucial. For low energy UHECR, the effects of intergalactic magnetic fields may be so strong 
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FIG. 1. Distribution of the arrival directions of 69 UHECR, represented by blue circles (o), with 
energy E > 5.5 x 10 19 eV reported by PAO in 2010, in the equatorial coordinates plotted using the 
Hammer projection. The solid red line represents the exposure boundary of the PAO experiment. 
The locations of Centaurus A (■) and M87 (♦) are also shown for reference. 

that the correlation of UHECR arrival directions with their sources can be completely erased. 
The energy cut 5.5 x 10 19 eV of the released data was chosen to be higher enough than the 
GZK cutoff so that their sources can be restricted within the GZK radius ~ 100 Mpc. This 
high value of the energy cut is good enough for correlation analysis. Thus, we use the full 
set of released PAO data for our analysis. 

From the Cen A dominance model, the expected arrival direction distribution can be 
obtained through the simulation. As an illustration, we show the arrival direction distribu- 
tions of 6,900 mock UHECR for two different values of the Cen A fraction, fc = 1.0 and 
fc = 0.3 with the same smearing angle 9 S = 30° in Fig. [2j From the observed UHECR 
set and the mock UHECR set, we obtain the observed and the expected CADDs/FEVDs 
as described in Eqs. (JBJ) and ((JJ), and calculate the KP statistic using Eq. (JSJ). To obtain 
the KP statistic with sufficient accuracy, we generate 10 6 mock UHECR events. Though we 
still have small fluctuations in the KP statistic for this number of mock events, the accuracy 
is enough for our purpose. Then, the probability is calculated by using Eq. (|JJ). We also 
checked the probability by direct Monte- Carlo simulation for several cases and confirmed 
that the Eq. (jUJ) is good enough. 

Fig. |3] shows CADDs and their CPDs of the PAO data, the isotropic distribution, and 
two cases of the Cen A dominance model. Close examination of CADD and its CPD is quite 
useful for understanding the results of statistical analysis and what causes the discrepancy 
between the data and the prediction of the model. For a single source, CADD is simply the 
distribution of angular distances of all UHECR from the source. The CADD of the isotropic 
distribution has a bell shape, reflecting the relative location of Cen A and the exposure 
function of the PAO experiment. On the contrary, the CADD of the PAO data has a 
distinguished peak at small angles, which means that there is an excess of observed UHECR 
near Cen A compared to the isotropic distribution. The KP test on CADD indicates that 
the probability that this excess of the PAO data is obtained from the isotropic distribution 
by chance is 1.3 x 10 -3 . Thus, this excess can be attributed to the Cen A contribution. 

The Cen A dominance model has two parameters, the Cen A fraction fc and the smearing 
angle 9 S , which can be used to fit this excess peak and the remaining bell-shaped body. The 
Cen A fraction fc affects the height of the peak at small angles relative to the bell-shaped 
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FIG. 2. Distributions of the mock UHECR arrival directions (6900 events, represented by small 
red dots) for the PAO experiment, obtained from the Cen A dominance model for two different 
values of the Cen A fraction fc = 1.0 and fc = 0.3 with the smearing angle 9 S = 30°. 




Angular Distance from Cen A Angular Distance from Cen A 



FIG. 3. CADDs (left panel) and their CPDs (right panel) of the PAO data, the isotropic distribu- 
tion, and the Cen A dominance model with the parameter sets (fc = 0.1, 6 S = 10°) and (fc = 0.3, 
9 S = 30°). Vertical lines in the right panel represent the sizes and the locations of the maximum 
differences between the CPD of the PAO data and those of the models considered. 
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FIG. 4. Cen A fraction (fc) and smearing angle (9 S ) dependence of probabilities that CADD (left 
panel) and FEVD (right panel) of the PAO data are obtained from the Cen A dominance model. 
The white + marks the parameter values of the maximum probability. 

body at larger angles. The smearing angle 8 S changes the position and width of the peak. 
As 9 S increases, the Cen A contribution tends to be isotropic and the peak merges into 
the bell-shaped body. Fig. H] shows the probability that CADD or FEVD of the observed 
PAO data is obtained from the Cen A dominance model, as a function of two parameters 
fc and 9 S . Compared to the CADD method, the FEVD method seems to becomes rather 
insensitive for small Cen A fraction or large smearing angles where the distribution tends 
to be isotropic. At angles smaller than 30°, both the CADD and the FEVD methods give 
the similar and consistent results. Thus, we mainly use the results of the CADD method 
to extract our conclusions. For the CADD method, the best fit parameters at which the 
probability reaches its maximum value are fc ~ 0.1 and 9 S = 5°, giving the maximum 
probability P max = 0.29. We don't have the la band (P > 0.32) in the whole parameter 
ranges. However, this is not so disappointing, considering the simplicity of the Cen A 
dominance model that the distribution of UHECR which are not contributed by Cen A is 
assumed to be isotropic. As you can notice from the CPDs in Fig. [21 the parameters in the 
2cr band fit the CPD at small angles well, but fitting gets poor at large angles where the 
isotropic contribution dominates. This makes the over-all fitting a little bit poor. The main 
cause of this poor fit is the big dip observed at angles 30°-40°, as seen in the CADD of 
the PAO data in Fig. [3j The big dip means the void region of UHECR at angular distance 
30°-40° from Cen A, which is actually observed at the lower right region near the south 
pole in the skymap of Fig. [TJ This departure from isotropy due to the void region makes 
the assumption of isotropic background in the Cen A dominance model a little bit poor 
one, though it is still a valid first approximation. Actually, this makes the best fit value 
of the smearing angle smaller than the naively (or more appropriately) expected value 10°~ 
20° at which the peak of CADD is observed. Considering this uncertainty concerning the 
assumption of isotropic background, we extend the viable parameter range for the Cen A 
dominance model to the parameter range of la band [P > 0.046): 0.05 < fc < 0.15 and 
0° ^ 9 S < 20°. In Section [Til we mentioned that fc is different from /c,pao, the fraction 
of Cen A contribution as observed at the PAO experiment. Because Cen A is located at 
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the central region of the PAO exposure, /c,pao is large than f c . /c.pao can be obtained 
as a function of fc and 8 S when we do the Monte-Carlo simulation. In terms of /c,pao> 
the best-fit value is /c,pao = 0.15 and the 2a band is 0.09 < /c,pao ^ 0.25, which means 
that among 69 UHECR observed by PAO, about 10 (2a band is 6 ~ 17) UHECR can be 
attributed to Cen A contribution. 



V. DISCUSSION AND CONCLUSION 

Let us discuss the implications of the results obtained in the previous section. First of 
all, though we started with the hypothesis that Cen A is a dominant source of UHECR, 
what we can actually prove through our statistical analysis is that there is a strong source 
of UHECR in the direction of Cen A. To confirm that the actual source is Cen A, other 
evidences, e.g., from the acceleration mechanism, energetics, or the energy spectrum, are 
needed. 

As seen in the previous section, the major source of uncertainty in our statistical analysis 
is the assumption that the background contribution is isotropic. We used the UHECR data 
with energies larger than 5.5 x 10 19 eV, which is larger than the GZK cutoff. Thus, their 
sources are believed to be mostly located within the GZK radius ~ 100 Mpc, where the 
matter distribution is not isotropic. One important issue in the UHECR arrival directions 
is whether the UHECR sources trace the matter distribution. In this regard, it is important 
to check the existence of the correlation between the UHECR arrival directions and the 
matter distribution within the GZK radius. There have been several studies on this, and 
the existence of correlation is not yet conclusive (8-10]. Thus, we stick to the isotropic 



background instead of modeling the matter-tracing background. Concerning the background 
contribution, we again draw your attention to the fact that the CADD of the PAO data 
depicted in Fig. [3] shows the large deficit at 30°-40° bin and another excess at 100°-130° 
range compared to the isotropic distribution. The deficit is due to the void region near 
the south pole in the PAO data. The chance probability that this void region is obtained 
from the isotropic distribution is about 5.0 x 10 -3 by the same CADD method and KP 
test. Not only the excess in the Cen A direction but also this void region makes the PAO 
data anisotropic. Another excess may be due to another point sources. We have examined 
the possibility of the existence of another point source at the region of angular distance 
100°-130° from Cen A, but found that this excess is consistent with isotropy. 

The importance of searching for the point sources of UHECR is that it is the starting 
point of cosmic ray astronomy. A known point source of UHECR can be used to probe the 
intergalactic magnetic fields in the vicinity of the source. For Cen A, this kind of study was 



done in Ref . [26[ . Here, we provide the order estimate from our results for comparison. Once 
we accept the results in the previous section that the fraction /c,pao = 0.15 of the observed 
UHECR by PAO is contributed by Cen A and their average deflection angle is around 10°, 
we can estimate the rms magnitude of the magnetic fields in the vicinity of Cen A. If we 
assume that the intergalactic magnetic field is composed of random patches of magnetic 



fields, the rough estimate for the deflection angle is given by [1| 



io 20 evy v iom p c V V 10 " 9G 

where E is the energy of the cosmic ray particle, d is the size of the magnetic field extension, 
£ c is the average size of patches, and B is the magnetic field strength. If we select 10 nearest 
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UHECR from Cen A, the average energy is E ~ 7.0 x 10 19 eV, and the average defection 
angle is 68 ~ 10° (The maximum deflection angle is 66 ~ 15°.). Inserting these values, Z = 1 
(the proton), and the distance d ~ 4Mpc, we obtain an estimate B ~ 14 (^/Mpc)" 1 / 2 nG. 
This gives the rough strength of the intergalactic magnetic field in the vicinity of Cen A. 

In conclusion, we examined the possibility that Cen A is a dominant source of UHECR 
observed by PAO, by using the statistical analysis of the arrival direction distribution. We 
set up the Cen A dominance model for the UHECR sources, in which Cen A contributes 
the fraction fc of the whole UHECR with energy above 5.5 x 10 19 eV and the isotropic 
background contributes the remaining 1—fc fraction. The effect of the intergalactic magnetic 
fields on the bending of the trajectory of UHECR originated from Cen A is parameterized by 
the gaussian smearing angle 8 S . We adopted CADD and FEVD methods for the reduction of 
the arrival direction distribution, and the KP test to compare the observed and the expected 
CADD and FEVD and to estimate the significance level of the similarity. We observed the 
excess of UHECR in the Cen A direction in CADD. Then we tried to fit the CADD of 
the PAO data by varying two parameters fc and 9 S of the Cen A dominance model. The 
best-fit parameter values are fc ~ 0.1 (/c,pao ~ 0.15) and 9 S = 5° with the probability 
P = 0.29. Considering the uncertainty concerning the assumption of isotropic background 
in the Cen A dominance model, we extend the viable parameter ranges to the 2a band, 
0-09 ^ /c,pao ^ 0.25 and 0° < 9 S < 20°. It supports the existence of a point source in 
the direction of Cen A, which is extended by the action of intergalactic magnetic fields. If 
Cen A is actually the source responsible for the observed excess, the intergalactic magnetic 
field in the vicinity of Cen A is estimated to be B ~ 14 (£ C /Mpc) _1//2 nG from the average 
deflection angle of the excess UHECR. 
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